Biomedical Letters ISSN 2410-955x- An International Biannually 


Research article i 
ages 25-30 


2023 | Volume 9 | iSSi 
b d 


ARTICLEINFO 


Received 
January 04, 2023 
Revised 

February 25, 2023 
Accepted 

March 27, 2023 


*Corresponding Author 
Samar Shahzadi 


E-mail 
samar.iiui@gmail.com 


Keywords 
Bioinformatics 


In silico structure 
prediction 


Hepatitis B virus 
Genotype A2 
HBeAg 


How to Cite 


Shahzadi S, Tahir Q, Khan 
MI, Raza A, Khan MF, 
Larra SA. In silico 
structure prediction and 
molecular docking 
analyses to reveal 
potential binding domain 
of Hepatitis B virus 
genotype A2. Biomedical 
Letters 2023; 9(1):25-30. 


In silico struct diction and 
molecular do alyses to 
reveal potential binding domain 
of Hepatitis B virus genotype A2 


Samar Shahzadit”, Qanta Tahir'?, Muhammad Imran 
Khan?, Ali Raza?, Muhammad Fardeen Khan’, Sajjad 
Ahmad Larra’* 


1Department of Bioinformatics, Institute of Biochemistry, Biotechnology 
and Bioinformatics, The Islamia University of Bahawalpur, Bahawalpur, 
Pakistan 
?Department of Microbiology and Molecular Genetics, University of 
Okara, Okara, Pakistan 

i y, Institute of Biochemistry, Biotechnology 
and Bioinformatics, The Islamia University of Bahawalpur, Bahawalpur, 
Pakistan 
4 National Institute for Biotechnology and Genetic Engineering (NIBGE), 
Faisalabad, Pakistan 


Abstract 

Hepatitis B Virus (HBV) infects the hepatocytes to cause 
serious liver diseases. HBeAg regulates the response of 
immune system to the intracellular capsid act as T-cell 
_tolerogen.The immune response regulation may predispose 
to chronicity during perinatal infections to prevent the 
severe liver injuries. Various in silico approaches including 
comparative modeling, threading approach and ab initio 
approach were employed for the prediction of 3D structures 
of the selected protein followed by the validation of the 
predicted structures through Errat, Procheck and Anolea. 
The predicted 3D structure of HBeAg revealed overall 
quality factor of 95.9184%. Interestingly, it was observed 
that only 1.97% residues were present in outlier region 
Oe ee a. While 98.03% in favored and allowed 
region. Molecular docking analyses 
were performed and the attempt was 
for the identification of novel ligands 
for HBeAg. The reported compound 
may regulate the activity and act as 
regulator of HBeAg. Interestingly, 
least binding energy of -7.1 Kcal/mol 
reported compound and high binding 


25 


Biomedical Letters 2023; 9(1):25-30 


affinity to predict the binding residues (Asp-51, Phe-53, Val-56, Arg-57, Met-95, Ala-98, 
Asn-103, Arg-111, Asp-112, Val-115, Val-118 and Asn-119). The function determination 
of the selected target protein is due to the identification of effective binding sites in 
protein structures. The reported compound may act as potent molecule and the 
predicted structure is reliable for the functional studies and structural insights. 


Scan QR code to 
Ir see this publication 

on your mobile 
Hepatitis B 
hepatocytes, 


infects 


(HBV) 
cell type of parenchymal 
tissues of liver to cause liver diseases. 
HBV life cycle is complex and has the 
ability to develop numerous antiviral 
agents. The life cycle of HBV has relaxed- 


virus 


circular partially double-stranded 
genomic DNA (rcDNA), converts into a 
molecular template DNA to enhance the 
viral RNA intermediate. Chronic infection 
results due to high stability of HBV which 
is poorly curable [1]. 250 million people 
are assumed to be chronically effected by 
HBV including one million deaths yearly. 
HBV uses numerous ways to tackle the 
host innate immunity to increase its 
replication includes taking advantage of 
the growing immune system of young 
children to smooth its persistence. It can 
also use maternal viral E antigen to guide 
the immunity of the offspring to support 
its perseverance after vertical 
transmission [2]. HBV is one of the 
familiar worldwide blood-borne pathogen. 
The chronic hepatitis B leads to an 
inactive carrier state which results in 
cirrhosis and fatal liver cancer. HBV 
surface-antigen vaccine is productive, but 
medicaments are now not curative [3]. 
HBeAG is an external core antigen 
protein with 214 amino acids. Main 
function of HBeAG is in regulating 
immune response to intracellular capsid. 
By having immune regulatory effect it 
acts as a T-cell tolerant to avert 
demolition of contaminated cells by 
cytotoxic T-cells. This immune regulation 
may incline to chronicity through 
perinatal infections and stop acute liver 
injury during adult infections [4]. 


This work is licensed under the Creative Commons Attribution 
Non-Commercial 4.0 International License. 


MbBeAg 1S a HBV replication Marker. 
HBeAg is the antigen of HBV and has the 
ability to cross the placenta towards 
specific insensitiveness of helper T cells 
to the capsid protein and HBeAg in 
newborns. HBeAg plays the role of 
tolerogen after birth as it is tolerated in 
utero. The HBeAg-positive mothers 
showed continual prenatal transmission 
while it is less incessant HBeAg-negative 


mothers. The genotypes and sub- 
genotypes of HBV may have perceptible 
geographical distribution leads _ to 


different mutations near to HBV genome 
coding for HBeAg. The distinct genotypes 
of HBV can be in charge of multiple 
infection natural history and different 
ways of transmission in children, present 
in numerous areas of the world, where 
different genotypes prevail [5]. 
Progressive improvements have been 
observed in immuno-informatics [6-10] 
and computational drug designing [11-21] 
from last decade. Numerous problems of 
biology have been resolved through 
applying various approaches of 
bioinformatics [14]. The present work 
demonstrates the molecular docking 
analyses to explore the compound against 
HBeAg. The ligands of vast structural 
entities and common structural features 
were explored. The experimental resolved 
3D structure of HBeAg through X-ray 
crystallography and Nuclear Magnetic 
Resonance (NMR) was not available yet. 
The 3D reliable structure of HBeAg was 
modelled by applying the crystal 
structure. 


Materials and Methods 
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The canonical sequence of HBeAg having 
accession number 091532 was retrieved 
from Uniprot Knowledgebase database in 
FASTA format. In present effort, the 3D 
structure prediction of the target protein 
was performed followed by the molecular 
docking studies on DELL Core-I-7 
workstation. The amino acid sequences of 
the target protein was retrieved and 
subjected to BLASTp to identify the 
suitable template by utilizing the Protein 
Data Bank (PDB) [22]. MODELLER 9.15 
[23], an automated program of protein 
modeling was used for the three 
dimensional (3D) structure prediction of 
the target protein by satisfying the spatial 
restraints. The retrieved templates were 
utilized to predict the structures through 
homology modeling approach [24, 25]. 
Various evaluation tools such as ERRAT 
[26], ProCheck [27], Rampage [28] and 
Anolea [29] were employed to evaluate 
the generated models. Moreover, the 
predicted structures of the target protein 
were further evaluated by using 
MolProbity evaluation tool [30]. 

ZINC library was utilized to analyze and 
evaluate the binding pockets of the target 
protein. The ligands were used and 
protein-ligand molecular docking 
analyses were performed by utilizing 
AutoDock tools [31]. The drug like 
properties including number of rotatable 
bonds, H-bond acceptors and H-bond 
donors were calculated by utilizing 
PubChem [32]. The Lipinski’s rule of five 
[33] was calculated for the selected 
ligand by utilizing mCule servers. The 
mutagenesis and carcinogenicity of the 
selected compounds were also calculated 
and no carcinogenic and mutagenic risks 
were observed. The objective to perform 
the molecular docking analyses was to 
identify the binding pattern of the target 
protein against the selected ligands. The 
geometrical optimization of the predicted 
3D structure of the target protein and the 
selected ligands were performed by using 
UCSC Chimera and ChemDraw Ultra 
respectively. The generated docked 
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complexes were visualized and analyzed 
by employing the VMD, PyMol, Ligplot 
and Chimera 1.6. The absorption, 
distribution, metabolism, excretion, and 
toxicity (ADMET) properties of the 
selected ligands were analyzed by 
employing the admetSAR. 


Results and Discussion 


The aim of this work was depends on the 
relation of Hepatitis B virus and in silico 
analyses to identify the binding pockets. 
The aligned templates against the 
selected target protein having E-value, 
maximum score, identity, and query 
coverage were selected for homology 
modeling (Table 1). The scrutinized 
suitable templates were employed for 3D 
structure prediction of the selected 
protein. The scrutinized templates 
showed reliable evaluation analyses. The 
3D predicted structures were modeled 
comparatively by utilizing the crystal 
structures of the selected templates 
retrieved from PDB. 

The utilized evaluation tools revealed the 
reliability and efficacy of the selected 3D 
predicted structures (Fig. 1 A-C). The 
favored, allowed and outlier regions of 
the predicted structures were calculated 
in Ramachandran plot. 98.03% residues 
were observed in favored and allowed 
region however, only 1.97% residues 
were observed in outlier region. The 
overall quality factor was observed 
95.91% for the predicted structure and 
seems reliable for further analyses. 
Moreover, the molecular docking studies 
of the reported ligands revealed 
fluctuation and variation in the binding 
energy. Initially, the molecular docking 
analyses were done having twenty poses, 
and 90 runs were saved, out of which the 
suitable effective poses of the complexes 
having least binding energy were 
selected. It was observed that the 
reported ligand effectively binds to the 
selected target protein (Table 2). 
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Table 1: The selected templates for External Core Antigen sorted by their E-values, overall 


quality, max score, query coverage and identity 


Description Maximum Total score Query Per E-value 
score coverage identity 

Chain B capsid 353 353 83% 94.94% 3e-125 
protein 
Hepatitis e-antigen 313 313 74% 94.34% 3e-110 
F97L Hepatitis B 302 302 69% 97.32% 2e-105 
protein 
CryoEM 310 310 69% 100% 2e-108 
Hepatitis B viral 307 307 69% 100% 7e-108 
capsid 


The molecular docking analyses were 
analyzed on the basis binding affinity, 
properties of the selected ligand, drug 
properties and least binding energy 
(Table 3). The selected ligand was 
observed as a cyclic compound (Fig. 2) 
having important and acceptable 
biological properties. The reported ligand 
may be considered as potent anti HBV 
agents by targeting HBeAg. 


Table 2: Molecular docking analyses 
including ligand efficiency and binding 
residues of the selected target protein 


Properties HBeAg 
Final intermolecular energy -10.87 
(kcal/mol) 
Estimated free energy of -7.1 
binding (kcal/mol) 
Estimated inhibition constant, Ki 37.76 
(uM) 
Unbound system’s energy -0.52 
(kcal/mol) 
Torsional free energy (kcal/mol) 4.04 
Ligand efficiency -0.47 
Binding residues Asp-51, Phe- 
53, Val-56, 
Arg-57, Met- 
95, Ala-98, 
Asn-103, Arg- 
111, Asp-112, 
Val-115, Val- 
118 and Asn- 
119 


The molecular docking tool was employed 
and top ranked complexes of the selected 
target protein of HBV with least binding 
energies were chosen for further 
analyses. Different variation in least 
binding energy was observed however, 
the stability of the selected ligand 
depends on the conserved binding 
affinities of the docked complexes. 
Interestingly, it was observed that the 
docked complexes of the selected target 
protein with the selected ligand showed 
reliable results by satisfying the selected 
parameters and filters. It was further 
revealed that the reported ligand 
molecule binds at the conserved region of 
the selected target protein leads to 
explore the binding residues. Extensive in 
silico analyses suggested that the 
observed binding residues and their 
combination may lead to the least binding 
energy of the docked complexes of the 
selected protein and reported ligand 
molecules (Fig. 3). 

The drug design process is a costly and 
time consuming [15]. Therefore, various 
computational approaches and 
techniques have been applied to design 
compounds [17]. The emergence of 
bioinformatics approaches have 


significance in decreasing the required 
time with minimum side effects [34]. 
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Fig. 1: 3D structure of the target protein HBeAg A) Ribbon structure B) Surface view C) Surface 


view having 80% transparency. 


The structures of the selected ligands 
were evaluated for their oral 
bioavailability and efficacy [35]. ADMET 
properties of the selected ligand were 
analyzed. Various mathematical models 
including carcinogens, acute oral toxicity, 
Ames toxicity, cytochrome P450 2D6 
inhibition, honey bee toxicity, aqueous 
solubility [LogS], blood-brain barrier 
penetration, human intestinal absorption, 
fish toxicity and Caco2 permeability were 


calculated. Numerous toxicities were 
predicted (Table 3). The observed 
toxicities help to evaluate the 


metabolites, intermediates and pollutants 
[12]. 


Table 3: The drug properties of the selected 
ligand (MLCB-0696) 


Ligand properties MLCB- 
0696 
Hydrogen bond donor 08 
cLogP -1.16 
Hydrogen bond acceptor 10 
Molecular weight (g/mol) 515.64 
Rotatable bonds 03 
Blood-brain barrier (BBB) 0.9000 
(probability) 
Human intestinal absorption 0.6995 
(HIA) (probability) 
Caco2 permeability (probability) 0.8570 
CYP450 2D6 inhibitor 0.9366 
(probability) 
Carcinogens (probability) 0.9400 
Acute oral toxicity (probability) 0.4931 
Aqueous solubility (LogS) -2.935 
Fish toxicity (LC50, mg/L) 0.5828 
Honey bee toxicity (HBT) 0.7079 


(probability) 


The predicted aqueous solubility (at 
25°C) of the selected ligand showed 
solubility of the compound in the water. 
The reported compound showed less 
LogP value leads to follows the Lipinki’s 
rule of five. It was observed that the 
compound could be easily be absorbed by 
the human intestine evaluated through 
mathematical model of intestinal human 
absorption. Toxicity and carcinogenicity 
risk assessment were also calculated, and 
it was observed that the selected ligand 
was non-carcinogenic. 


Extensive literature survey and in silico 
analyses suggested that the suitable 
compound must be the one 


Fig. 2: The 2D structure of the selected ligand 
(MLCB-0696) 


that undergoes from the selected 
parameters to satisfy the drug properties, 
least binding energy and effective binding 
affinity values. By applying the selected 
parameters, it is suggested that the 
reported compound has potential to use 
against HBV by targeting HBeAg. The 
generated molecular docking results 
suggested that the observed binding 
residues (Asp-51, Phe-53, Val-56, Arg-57, 
Met-95, Ala-98, Asn-103, Arg-111, Asp- 
112, Val-115, Val-118 and Asn-119) were 
crucial. 


Conclusion 


In conclusion, the reported ligand showed 
efficacy against HBV by targeting HBeAg 
through computational analyses. 
Extensive in silico analyses of HBeAg 
showed higher efficacy and probability 
based on used parameters and least 
binding energy. The potential interacting 
residues (Asp-51, Phe-53, Val-56, Arg-57, 
Met-95, Ala-98, Asn-103, Arg-111, Asp- 
112, Val-115, Val-118 and Asn-119) 
identified by molecular docking analyses 
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may be significant for site-directed 
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